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We describe the quantum mechanical spreading of a Gaussian wave packet by means of the 
semiclassical WKB approximation of Berry and Balazs. We find that the time scale r on which 
this approximation breaks down in a chaotic system is larger than the Ehrenfest times considered 
previously. In one dimension r — ^X~^ ln{A/h), with A the Lyapunov exponent and A a typical 
classical action. 
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According to Ehrenfest 's theorem the propaga- 
tion of a quantum mechanical wave packet is described 
for short times by classical equations of motion. The 
time scale at which this correspondence between quan- 
tum and classical dynamics breaks down is called the 
Ehrenfest time. If the classical dynamics is chaotic with 
Lyapunov exponent A, then the Ehrenfest time r is of 
order A~^ln(^/?i) (with A a typical classical action of 
the dynamical system) Q . There is actually more than a 
single Ehrenfest time, corresponding to different types of 
semiclassical approximations. Although they differ only 
by a numerical coefficient, Ti = CiX~^ ln{A/h), the struc- 
ture of the wave function changes qualitatively from one 
time scale to the next. 

Up to a time ti, with ci — 1/6, the initial coherent 
state will retain its Gaussian form with vanishing error 
in the limit h —> . For longer times up to T2 , with 
C2 = 1/2, the uncertainty in position and momentum of 
the particle remains small but the phase space structure 
of the wave packet deviates strongly from a Gaussian. For 
times greater than T2 the wave function no longer has the 
form of a wave packet (this is the "mixing regime" of Rcfs. 
lH,|j), but up to a time ra it can still be described semi- 
classically by the time-dependent WKB approximation of 
Berry and Balazs j7|. As we will show in this paper, the 
WKB representation implies C3 — 7/6 for a single degree 
of freedom (with simple generalizations for higher dimen- 
sions) . This is larger than the value C3 = 2/3 obtained by 
Bouzouina and Robert from a different semiclassical 
approximation. 

Let us start with the Gaussian one-dimensional wave 
packet 
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Initially f3{t = 0) = and a{t = 0) = pf/L, where pp 
and L are the typical classical momentum and length. 
The typical classical action is ^ = ppL. The parame- 
ters xo{t),po{t) follow the classical trajectory for h A. 
We will measure momentum and coordinate in units of 
Pf and L, respectively, so that a(0) = 1 and A = 1. 



For chaotic dynamics with Lyapunov exponent A one has 
a{t) oc exp(— 2At), hence a <C 1 for t ^ 1/A. 

To describe the time evolution in phase space we con- 
sider the Wigner function 
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The wave packet is centered at xo{t),po{t) and for 
a{t) <C 1 becomes highly elongated and tilted with slope 
Ap/Ax « p. It has length = y/K(TT]P)/a and width 
lj_ — y/haJJT+fW), so that the area in phase space is 
conserved exactly, l\\l± = h. The Gaussian quantum 
wave packet satisfies the classical Liouville theorem. 

The Gaussian form (|^) takes into account the elonga- 
tion of the wave packet, but not the curvature that de- 
velops in time and results in a bending of the packet. To 
describe the curvature we add an imaginary cubic term 
in the exponent in Eq. (|^), 
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(For simplicity we have put xq = 0.) The cubic term 
leads to an appreciable phase shift over a length l\\ ~ 
(hjaYI"^ when /h){h/ aY^"^ > 1, hence when a(t) < 

For a < ^Yi^ 

gncr function takes again a 
simple form, in terms of the Airy function Ai: 
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One can check that W{x,p) — > 5{x)5{p — pi^) when h Q 
(at fixed a), by means of the identity lim^^o Ai(z/£)/e — 
Y^(5(z). At finite h the wave packet is extended along 
the curved line p = pq + l3x-\-^x'^ . Since p,po,x are of or- 
der unity, the two parameters /3 and 7 are of order unity 
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as well (in contrast to a, which is ^ 1). The transverse 
width is of order 
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The length of the packet remains at Z|| « y/h{l + fW)/a. 
Since now 3> h, the Liouville theorem no longer 
holds. 

To obtain the Ehrenfest time we parameterize time as 
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The classical limit for a chaotic system means h ^ 0, 
t ^ oo at fixed c. Different coefficients c follow from dif- 
ferent semiclassical approximations. If we use the Gaus- 
sian wave packet (|l|), without the cubic term to account 
for the curvature, then we need a{t) ^ h^^'^y^l'^ . Since 
a oc e"^"^* (X Ti''^ we need c < 1/6. The upper limit of c 
gives the first Ehrenfest time t\ — ^A"^ ln(l/?i). 

The classical limit can be reached for longer times if we 
use the wave packet (||), including the cubic term. The 
dimensions of the packet for t>T\ scale with fi as 
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For c < 1/2 the length of the packet approaches zero in 
the classical limit. This upper limit of c gives the second 
Ehrenfest time T2 = ^A"^ ln(l/?i). 

For t > T2 the length of the wave packet exceeds the 
size of the system and is no longer small compared to the 
radius of curvature. For these large times we may adopt 
the semiclassical WKB approximation of Berry and Bal- 
azs [0. Consider a curve in phase space and a phase 
space distribution p{p{x),x). Both p and p evolve in ac- 
cordance with classical equations of motion. For t > T2 
the function p{x) is multivalued with an exponentially 
large number of branches ~ exp[A(t — T2)]. The quantum 
wave function in this "mixing" regime has the form 
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The summation over k accounts for the different branches 
of the multivalued function p{x). The two functions / 
and a are related for ?i — s- to p and p by the correspon- 
dence principle, 



^ = p{.x) , / = \/p{p7x) ■ 
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An explicit description of the evolution of the wave func- 
tion (|^) for quantum maps may be found in Rcf. 

Near the point Xh at which p{x) bifurcates into two 
branches, one has p = pt ± a^/x — Xb, p — b/ ^Jx — Xh- 
The wave function there is 

= {hjaf^ b^'^kli (^{a/hf/^ (x - Xb)') e*^"'^ (10) 



up to an overall phase. The phase difference between the 
bifurcation points can be determined from Eqs. and 
(^). Because the curve p{x) is not closed, there is no 
analog of the Bohr-Sommerfcld quantization rule. 

The Wigner function corresponding to the wave func- 
tion (||), being quadratic in 5*, contains both diago- 
nal {Wkk oc l/fcp) and oscillating nondiagonal {Wkm oc 
flfrn) contributions. Far from bifurcations the diagonal 
contributions to the Wigner function read 
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We have made a Taylor expansion of (t(x ± y/2) and ne- 
glected the difference between f{x ± y/2) and f{x). 

If we parameterize time as in Eq. (^ we have for both 
Z|| and the same scaling with U as in Eq. (^. The range 
of validity of Eq. (||) is limited by the condition that dif- 
ferent branches should be distinguishable. This requires 
that the different parts of the curve p{x) in phase space 
should not get closer than l^. Their spacing is of order 
1/Z|| (assuming a uniform filling of phase space), hence 
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The upper limit of 7/6 for c leads to the third Ehrenfest 
time, 
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The third derivative a'" in Eq. ( |11| ) vanishes at the 
points of inflection of the curve p{x). In order to find 
the Wigner function there one should expand (7(2; ± y/2) 
up to terms of order y^ . This leads to a different scaling 
li_ oc fi^^'^ of the width of the Wigner function near the 
inflection points. Because these are isolated points, they 
will not contribute to the matrix elements of nonsingu- 
lar operators (containing only smooth functions of x and 
p). This different scaling should therefore not affect the 
Ehrenfest time (|l^). 

The nondiagonal contributions Wkm to the Wigner 
function lead to the "ghost curves" discussed in Ref. 
The Wigner function near these curves is given by the 
same Airy function as in Eq. (ITTl), but in addition ac- 
quires a strongly oscillating factor. Due to these oscil- 
lations the nondiagonal terms do not contribute to the 
matrix elements of nonsingular operators. (They may 
play a role in the decoherence by the environment Ipof.) 
At t > T3 the ghost curves merge with the (multivalued) 
curve p{x) and become indistinguishable. 

The time scale (13) for the breakdown of the 
WKB approximation is greater than the Ehrenfest time 
|A~^ln(l/?i) in the mixed regime obtained in Ref. 
That shorter time scale may signal the breakdown of the 



2 



series expansion ak{x) —>■ J2j=o^kjix)h^ ■ Then Eq. (||) 
would no longer hold, but for t < T3 the representa- 
tion (H) with a renormahzed function ak{x) would still 
be valid. 

So far we have discussed a one-dimensional (ID) 
chaotic system, which in general can be represented by an 
area preserving map . A familiar example is the kicked 
rotator pj| ]. For mesoscopic quantum dots, however, 
a more relevant model is the d-dimensional {d = 2, 3) 
Schrodinger equation with a smooth potential V{'r). The 
Gaussian wave packet then takes the form 
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Here S is the action for the classical trajectory ro(t) 
and we have defined pb ~ niVQ, x = f — ra,Qn = 
(3in + iain- As before, we rescale momentum and co- 
ordinate such that the typical classical action A — 1. 
Initially, Qn — iSin- Similarly to the one-dimensional 
case, ain defines the form of the packet in coordinate 
space and = Api/Ax„ give the angles in phase space. 
Substituting the wave function (|lj) into the Schrodinger 
equation one finds Newton's equation of motion for vq. 
The spreading of the wave packet in phase space is de- 
scribed by 
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This is the equation describing the spreading in phase 
space of a small Gaussian bunch of classical particles. 

The Wigner function corresponding to the wave 
function (|lj) has the Gaussian form W oc 
exp{-QiMinQn/h), where Q = (r-ro,p-po) is a 
vector in 2d dimensional phase space. The d Lyapunov 
exponents Xi (i = 1,2, ...d) govern the large-time behav- 
ior of the eigenvalues rrii = l/m2d-i+i oc exp(2Ait) of the 
real symmetric matrix M. Because of energy conserva- 
tion one Lyapunov exponent vanishes. We order the A's 
from large to small, so that Ai is the largest and A^ — 0. 

The wave packet remains Gaussian (preserving the vol- 
ume cx h"^ in phase space) until the curvature starts to 
play a role (via a cubic term in the action). The corre- 
sponding Ehrenfest time ti = ^Aj~^ ln(l/fi.) is the same 
as in ID, only now it is defined through the largest Lya- 
punov exponent Ai. The second Ehrenfest time, when 
the length of the packet exceeds the size of the system, 
also has the same form T2 = ^\^^ln{l/h). 

The third time T3 is different for d = 2, 3 from the ID 
case. Instead of Eq. (|7|) one now has 
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The longitudinal dimensions correspond to eigenval- 
ues rrii with 1 < i < d— 1, and the transverse dimensions 



l^^ to rrii with d + 2 < i < 2d. The two unit eigenvalues 
nid — rud+i — 1 contribute another factor ^/h each to 
the total volume V in phase space covered by the wave 
packet: 
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The available area Vmax is restricted to a shell of con- 
stant energy with thickness v^, hence Vmax oc \Jfi. We 
require V % Vmax for the semiclassical approximation, 
which leads to the Ehrenfest time 
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In conclusion, we examined different time scales = 
CiA"""^ ln(l/?i) for the breakdown of different types of 
semiclassical approximations. These Ehrenfest times dif- 
fer only by a numerical coefficient c^, which may seem 
insignificant. However, this difference is actually a signal 
of a different power law scaling with fi of the volume V 
in phase space covered by the wave packet. For short 
times Liouville's theorem dictates V cx ?i. For long times 
[parameterized as i = (c/A) ln(l/?i)] the WKB approxi- 
mation gives V cx fi 1^^'^ for a one-dimensional quantum 
map (such as the kicked rotator) and V cx f^J'^l^^'^l^^'^ 
for a d-dimensional conservative system. These different 
power laws reflect the fundamental change in the struc- 
ture of the wave function with increasing time and should 
therefore have observable consequences. Two possible 
applications are the Loschmidt echo |l^ and the quan- 
tum shot noise , where the Ehrenfest time plays a key 
role. 
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